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The multiconfiguration valence-bond 
method (VB) is applied to diatomic 
molecules using the Hartree-Fock (HF) 
atomic basis set. The hyperfine constant, 
Fermi contact term, is computed as a func- 
tion of the interatomic separation for the 
X^n ground state of ^^OH and X^X ground 
state of ^^^AgH"^ molecules. This study 
leads to a number of conclusions about the 
influence of correlation and polarization 
effects on the hyperfine structure of hydro- 
genic molecules. The calculated values of 
the Fermi contact term are found to agree 
within 1 % of the experimental values 
wherever available. 
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1. Introduction 



In this work, the hyperfine Fermi contact term as a 
function of internuclear distance R is computed for the 
X'n state of OH and the X'2 state of AgH^ The all 
electron ab initio valence-bond method [1] with 
Hartree-Fock (HF) and Sturm's [2] basis sets is used to 
calculate the electron spin density near the nuclei cor- 
rectly. This leads to an accurate determination of the 
molecular magnetic dipole coupling constants and in 
particular the molecular Fermi contact term, Ac. The 
method uses many-electron atomic wave functions to 
construct the molecular wave function and yields the 
correct asymptotic properties for the molecule, one of 
which is the convergence of the molecular hyperfine 
interaction parameters to the atomic values for large 
internuclear separations. 

Magnetic hyperfine parameters are very sensitive to 
the quality of the molecular wave function in general and 



to the spin polarization of atomic cores in particular. 
One of these parameters, the Fermi contact term, is 
proportional to the electron spin density at the position 
of the nuclei that have nonzero spin. When a molecule 
has electrons in open shells, it leads to one or more 
unpaired spins in valence orbitals and gives an un- 
balanced spin density at the nuclei. Computationally it 
implies different exchange potentials for electrons with 
spin up and spin down. The contact term, Ac, is pro- 
portional to the difference in the spin densities of 
electrons with the opposite direction of spin. The 
accurate determination of this parameter requires both 
correlation and polarization interactions in the model. 
The configuration interaction (CI) approach is used to 
treat these effects. We build the CI on the basis of the 
nonorthogonal many-electron atomic HF and Sturm's 
functions. 
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Hyperfine splittings have been observed in the optical 
spectra of diatomic molecules with very high accuracy 
Very precise theoretical calculations are required for 
comparison. One of these calculations was performed 
by Kristiansen and Veseth [3] for the OH molecule. A 
many-body perturbation theory was used to compute 
magnetic hyperfine parameters for the lowest vibra- 
tional levels of the ^11 ground state. Their disagreement 
with experiment does not exceed 2 % near the equi- 
librium distance /?e = 1.8342 a.u. (the atomic unit of 
length, a.u., is the Bohr radius ao and is approximately 
equal to 0.0529177 nm). 

In our calculations the Fermi contact term is deter- 
mined over a wide range of R (1.5 a.u. to 8 a.u.). The 
asymptotic value of Ac can be compared with our value 
at the largest internuclear separations. 

The only previous study of the molecular potentials 
of the AgH"^ molecule [4] was performed using a pseu- 
dopotential approach with its inherent integrated treat- 
ment of the core orbitals, and therefore does not give 
theoretical information about the AgH"^ hyperfine 
coupling constants. Our AdR) values were obtained for 
R between 2 a.u. and 1 8 a.u. and converge for large R to 
the atomic values. 



2. Theoretical and Computational Details 

The total electronic wave function ^ab of the 
molecule is introduced as a linear combination of 
Slater's determinants det^, corresponding to various 
configurations a of atoms A and B. That is 



^AB = ZjCadeta 



(1) 



where the Ca are obtained by solving a generalized 
eigenvalue matrix problem, described by the equation 



^AB C — SabC , 



(2) 



where Hab is the Hamiltonian of the molecule AB. The 
right-hand-side of Eq. (2) includes the nonorthogonality 
matrix S'ab, which describes an overlap between deter- 
minants <detjdet^> and it is given as 



(^AB)a3 = <detjdet^> = (DaaD^^) ^^Dap, 



(3) 



particular one-electron orbitals a^ and (3^, the overlap 
matrix elements 5"/ are used to describe the one-elec- 
tron density matrix p?'^ (x,x') of the molecule as 

p?'^(x,x') = (Z),,Z)33)-''^Z),3 

N 

^(s-')r/-Mx)-^j'(.x'), (4) 

where the <^(x) are the one-electron wave functions and 
X denotes both coordinates and spin of the electron. 
Moreover, the two-electron density matrix is 

P?'^(xi,X2lx{,xD = (D^aD^^)-''^ •D.p 

N N 

E E£>s!l/ • <^. (xi)<^;(xo<^. fe)<^;(x3, (s) 

ij k,l 

where 

D^h., =£>„p • [(5-'),^f* • {S-')i:hS-')^i^ ■ (5-')£f*] 

(6) 

Finally, the Fermi contact term can be expressed in 
terms of the total spin densities p"'^(T) and p"'^(nI) at 
each nuclear site for the electrons with spin pointed up 
(T) and down (i). The CI form of Ac is introduced as 



8tt 
3 



A. = 2lXo gi Miv E ^c„C3[p"'^(T) - p"'^(i)] (7) 



D^^ = det\<oii\^i > ... <a/ip^- > ... <aiviPiv> I, and N 
is the total number of electrons in the molecule. For the 



where gi denotes the nuclear g -factor, ^lcn is the nuclear 
magneton, and ^Xq is the Bohr magneton. 

In our model the CI expansion for OH includes con- 
figurations obtained from the occupied \s^, 2s^, and Ip"^ 
shells of oxygen and the 1^ shell of hydrogen, as well as 
configurations created by excitations of these electrons 
into the 3^, 3p, 3d, and 4^ states for O and the 2s, 2p, 
3^,3/7, and 3d for H. The total number of configurations 
was equal to 238. To describe the virtual excited states 
we use localized Sturm's functions. The result of our 
computation of the A^iR) for the oxygen and hydrogen 
atom in the ^^OH molecule is shown in Fig. 1. The 
agreement of the calculated values over the wide range 
of internuclear separations R (1.5 a.u. to 8 a.u.) is within 
less than 1 % of the precise experimental values, which 
are for the equilibrium distance /?e= 1.8342 a.u., 
Ac = -73.1258 MHz [5], and for the asymptotic values 
A(0) = -219.6 MHz and AfH) = 1420.40575 MHz 
[6, 7]. 
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Fig. 1. Fermi contact term Ac of the OH molecule as a function of 
the intemuclear separation (1 a.u. = 0.0529177 nm). 



Fig. 2. Fermi contact term Ac of the '°^AgH^ molecule as a function 
of the intemuclear separation (1 a.u. = 0.0529177 nm). 



The CI treatment of the AgH "^ molecule was based on 
the occupied 4p^4d^^ shells for the Ag"^ ion and the 
1^ shell for the H atom as well as the 5^, 5p, 5d, 6s 
and 2s, 2p, 3s, 3p, 3d virtual Sturm's orbitals for 
Ag"^ and H, respectively. Computed values of the Fermi 
contact term of ^^^AgH"^ molecule are presented in 
Fig. 2. The asymptotic values of Ac for H and Ag^ fit 
with about 1 % disagreement to the well known experi- 
mental value of the hyperfine coupling constant for the 
H atom, and to the expected zero value for the closed 
shell ion Ag^. 
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